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Abstract Water and water-mediated interactions determine thermodynamic 
and kinetics of protein folding, protein aggregation and self-assembly in con- 
fined spaces. To obtain insights into the role of water in the context of folding 
problems, we describe computer simulations of a few related model systems. 
The dynamics of collapse of eicosane shows that upon expulsion of water the 
linear hydrocarbon chain adopts an ordered helical hairpin structure with 
1.5 turns. The structure of dimer of eicosane molecules has two well ordered 
helical hairpins that are stacked perpendicular to each other. As a prelude 
to studying folding in confined spaces we used simulations to understand 
changes in hydrophobic and ionic interactions in nano droplets. Solvation of 
• hydrophobic and charged species change drastically in nano water droplets. 

Hydrophobic species are localized at the boundary. The tendency of ions to 
be at the boundary where water density is low increases as the charge den- 
sity decreases. Interaction between hydrophobic, polar, and charged residue 
are also profoundly altered in confined spaces. Using the results of computer 
00 ' simulations and accounting for loss of chain entropy upon confinement we 

argue and then demonstrate, using simulations in explicit water, that or- 
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dered states of generic amphophilic peptide sequences should be stabilized in 
cylindrical nanoporcs. 

Keywords confinement effects on protein stability • hydrophobic interac- 
tions • potentials of mean force ■ water in pores ■ collapse and association of 
hydrocarbon chains 



1 Introduction 

In the process of protein folding an ensemble of unfolded states reach compact 
folded structures. Protein folding has become a paradigm problem in molecu- 
lar biology, and has inspired a large number of studies using numerous meth- 
ods [551l42U47[H61l461l56j . The cited reviews have shown that a multi- faceted 
approach is needed to understand global aspects of protein folding. Although 
deciphering the global principles of folding cannot be achieved using brute 
computer simulations alone, specially designed computers that are capable 
of generating long trajectories, are starting to reveal molecular details of the 
folding of small single domain proteins and the role of water in well-defined 
systems (48] . Historically, considerable understanding of the nature of driv- 
ing force for protein self-assembly have been obtained simple model systems. 
In this review, we describe computer simulation results for simple systems, 
which illustrate the role water in protein collapse and folding in confined 
spaces. 

Typically hydrophobic residues are buried in the interior of folded pro- 
teins, and hence are shielded from water. The predominant driving force in 
protein folding is hydrophobic interaction between non-polar residues |54U15[ 
[3] , which is reasonable because globular proteins contain ~ 55% hydrophobic 
residues [7]. Thus, understanding the hydrophobic effect, and more generally 
water-mediated interactions, in the context of protein folding and aggregation 
[STllS^llfJS] holds the key in describing their self-assembly. Hydrophobic inter- 
action refers to the free energy increase upon transfer of non-polar species 
from non-polar solvents to water. The reluctance of water to allow rupture of 
the persistent but dynamically changing complex network prevents non-polar 
species from being easily solvated. To dissolve small (< 1 nm) non-polar so- 
lutes, a cavity of appropriate size and shape needs to be created in water, 
which is free energetically unfavorable. Water-mediated interactions between 
solute molecules play a crucial role in protein folding [15,32 . Interactions 
between small hydrophobic solutes like methane have served as useful stand- 
ins for those between amino acid side chains in aqueous solution 1 1 1 j . Such 
interactions, quantified using potentials of mean force (PMFs) averaged over 
solvent configurations, have been extensively studied in bulk water [431110) . 
However, globular proteins are polymeric and undergo collapse transition 
from unfolded states, and hence polymer chains that can adopt globular 
conformations in water are good model systems for probing the process of 
collapse. The first part of the perspective deals with new results for collapse 
of a single hydrocarbon chain and the association of two chains. These studies 
are models for describing the early stages of protein collapse and aggregation 
driven by hydrophobic forces as is the case for A/3-pcptides. 
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In the second part, we turn our attention to changes in molecular in- 
teractions between amino acid side chains in confined spaces. Besides their 
relevance in materials science, confinement effects play a major role in protein 
folding in vivo, e.g., in chaperonin-assistcd folding |57j . co-translational fold- 
ing in the exit tunnel of the ribosome [651169] . and dynamics in the crowded 
cellular environment {5511121128] . Several experimental [T71IT51I551I51 I5] . theoret- 
ical and computational studies 4,68, 25, 69, 13 , 33 , 67 have examined changes 
in the confinement-induced stability of proteins. General theoretical consider- 
ations and explicit simulations using coarse-grained models for polypeptide 
chains have shown that the subtle balance of many factors, including hy- 
drophobic and ionic interactions in confined water, entropic restrictions of 
the conformations of polypeptide chains, and specific interactions between 
amino acid residues and the confining boundaries determine the stability of 
confined proteins [41 IT3"1I4*T] . In many cases, stabilization of the folded state 
of proteins in confined spaces relative to bulk solvent [TTlHBI fSllS] can be ex- 
plained by the entropic stabilization mechanism (ESM). According to ESM 
reduction in the allowed conformations of the unfolded states of proteins in 
confined leads to stabilization of the native state with respect to the bulk 
4 , 68 , 25 , 69 . However, confinement can also destabilize the folded state [TBI 
133] due to the alterations in hydrophobic interactions, which can result in a 
net attraction between the protein and the confining boundary [T7lfT8] . To 
understand confinement effects systematically we describe the changes in hy- 
drophobic and ionic interactions using model systems and side chains (SCs) 
of amino acid residues in spherical and cylindrical cavities, and expound on 
the implications of these findings for folding in confined spaces. 



2 Methods 

2.1 Hydrocarbon chains in water 

Models: We give a brief description of the simulations reported here for 
the dynamics of collapse and association of hydrocarbon chains. Intermolec- 
ular interactions are modeled using standard site-site potentials. The pair 
potentials include partial charges and Lennard- Jones (LJ) interactions. The 
charges are located on the centers of the atoms. Site-site interaction be- 
tween molecules i and j, with R denoting collectively the coordinates of the 
molecule, is given by 

atoms /~<\2 v^6 

V(R i ,R j )= Y ^L_^l + ^lL_ (i) 

where r a fj is the scalar distance between sites a and j3. With this standard 
parametrization the interaction between molecules is fully specified in terms 
of the LJ potentials and the partial charges q a . 

Water: We use the three-site rigid SPC/E model for water in which the 
hydrogen atoms are at a distance 0.1 nm from the oxygen site and the HOH 
angle is 109.47°. The charges on oxygen and hydrogen are qo = — 0.8476e 
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and qn = 0.4238e, respectively. The LJ parameters for the oxygen sites are 
to = 0.1548 kcal/mol and ao = 0.3166 nm. 

Hydrocarbon (HC): Following our previous studies, we used a united 
atom, site-site interaction model for the linear alkane chain CH% — (Ci?2) n ~ 
CH3 with n = 18. The intramolecular force field for the hydrocarbon chain 
can be decomposed into Vi ntr a = V2 + V3 + V& + V5. The first two terms 
correspond to harmonic bond stretching and harmonic bond angle poten- 
tial. The term V4 contains the torsion interactions involving four consecutive 
sites. The Vg term represents non-bonded interactions between sites that are 
separated by three or more sites along the chain. The non-bonded terms are 
modeled using Lennard- Jones potentials. The parameters for the HC force 
field are given elsewhere [50] . 

Simulation Details: The system consisted of 1995 (1942) solvent molecules 
plus 1 alkane chain molecule (2 alkane chain molecules). The equations of mo- 
tion were integrated using an iterated form of the Beeman algorithm with 
a time step of 1 fs. Periodic boundary conditions were imposed in all three 
spatial dimensions. The volume of the system was set so that the water den- 
sity was 1 g/cm 3 , and the temperature of the system was maintained at 25 C 
by Nose-Hoover thermostat. At the start of the simulations, the alkane chain 
was placed in contact with the solvent and the solvent molecules were allowed 
to adapt to the presence of the alkane before any motion of the alkane was 
permitted. In this way, the disruption of the solvent due to the presence of 
the solute was allowed to relax. Once the systems were stabilized, production 
runs were made for long enough duration so that both collapse of a single 
chain and association of two chains could be observed. 



2.2 Small molecules and amino acid side chains in water 

We used Metropolis Monte Carlo [37] simulations to study the energetics of 
aqueous solvation of methane molecules and the model ions M g+ and M g _ in 
spherical and cylindrical water-filled nanopores using the TIP3P model for 
water [35] and a unified atom representation for methane [21]. Model ions 
M q ± were created by adding charges of magnitude q + = \q~\ = 0.4e or e to 
the methane spheres, where e is the electronic charge [S3]. In order to assess 
the importance of charge density (Q) in affecting solvation we also studied, 
for the 3 nm spherical droplet, ions M* 9+ and M* g _ with q + = \q~ \ = e and 
twice the volume of the methane sphere. 

Simulations were performed at constant number of molecules N, total 
volume V and temperature T = 298 K, and thus sampled the canonical 
ensemble. The effective volume available to the water molecules is taken to 
be 

V cS = V - N S V S (2) 

where N s is the number of solute molecules in the pore, and V s is the ex- 
cluded volume due to each solute. The spherical droplets had diameters (D) 
in the range 1.0-4.0 nm. The cylindrical pores had diameters (D) in the 
range 0.8-20 nm, and lengths L varying from 1.2 to 8.3 nm. These volumes 
were bounded by hard walls, with the potential energy at the wall being 
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large enough (10 12 kJ/mol) to confine the system to the desired volume. To 
simulate confinement effects, no periodic boundary conditions are applied. 

In the spherical droplets, the free energy of a single solute molecule was 
calculated as a function of its radial distance from the center of the droplet 
using Monte Carlo simulations with umbrella sampling and the weighted 
histogram analysis method (WHAM) [29]. Umbrella sampling and WHAM 
were also used to calculate the PMFs between pairs of solutes in the spherical 
droplets. In the cylindrical pores, the same technique was used to calculate 
interaction free energies, — fc^T log P(r), where ks is Boltzmann's constant 
and P(r) is the probability of finding the two solutes a distance r apart. 
Because of the quasi one-dimensional nature of the confinement when r D, 
we do not subtract the free energy contribution — 2k bT log r which arises 
from the increase in phase space proportional to r 2 in spherically symmetric 
systems. Therefore, these profiles cannot be compared directly to PMFs in 
bulk or in the spherical droplets. 

Molecular dynamics simulations: The interaction between the side 
chain pairs ALA-PHE, SER-ASN and LYS-GLU was studied [ST] using molec- 
ular dynamics simulations with the CHARMM22 force held [5] and TIP3P 
water |22j . In each amino acid, only the side chain was retained. The carboxyl 
and amino groups were deleted and the C a atom was replaced by a hydrogen 
atom. Thus, for instance, the side chain of ALA was represented by methane 
and PHE by toluene. 

In the cylindrical pores, interaction free energies, -fc^T log P(r) (ks is 
Boltzmann's constant, r is the distance between the centers of mass of the 
SC pair, and P(r) is the probability of finding the two solutes at a separation 
r), were calculated at a fixed volume using using the adaptive biasing force 
(ABF) technique [131121] implemented in NAMD [HJ. Because of the quasi 
one-dimensional nature of the confinement when r D, we do not subtract 
the free energy contribution — 2fcsTlogr which arises from the increase in 
phase space proportional to r 2 in spherically symmetric systems. Therefore, 
these profiles cannot be directly compared to PMFs in bulk. 



3 Results and Discussion 

3.1 Hydrocarbon Chains in Water: 

Hydration of Octane: The probability of creating a large anisotropic cavity 
to accommodate a solute would decrease with the volume of the solute. Thus, 
one might expect that large solutes would distort the hydration structure at 
least locally. Such an expectation is borne out in the studies of water near 
a hydrophobic wall, which showed that close to the wall there is a loss of 
one hydrogen bond per water molecule 30,31 . A natural question is how 
the solvent responds to the presence of large apolar molecules with internal 
degrees of freedom like linear hydrocarbon chains CH3 — {CH%) n — CH3, 
which are interesting model systems for studying the early stages of collapse 
of proteins. As a result, a number of studies focusing on various aspects of 
hydrocarbon behavior in water have been reported [40ll9ll20] . We first discuss 



6 



the case of octane relegating the discussion of hydrocarbon collapse to the 
next section. 

Because linear hydrocarbon molecules are intrinsically flexible and can 
undergo rotational motion around the dihedral angles populating cis, trans, 
and gauche states, one may expect that a solute with the size of octane 
(n = 6) or dodecane (n = 10) could locally distort the structure of water. 
Indeed, if octane or dodecane were treated as a sphere the probability of 
creating a void large enough enough to accommodate them is exceedingly low. 
Molecular simulations studies of solubility of large hydrocarbon molecules 
revealed several interesting findings. (1) Using the SPC/E model for water 
and a coarse-grained hydrocarbon model described in [50j . Mountain and 
Thirumalai |39) showed that there are minimal changes in the tetrahcdral 
network of water structure even near the hydration layer of octane. Although 
changes in water structure was not reported in simulations of dodecane [62] it 
was found that the extended structure was favored implying that the network 
of hydrogen bonds is sufficiently plastic to accommodate a large solute as 
large as dodecane. (2) Comparison of the distributions, P(cos8), {9 is the 
angle between three near neighbor oxygen atoms) between water in the bulk 
and oxygen atoms within 0.5 nm of at least one of the carbon atoms in octane 
shows peaks around the tetrahedral value (cos8 ~ 5) [35], which shows that 
a cavity of an appropriate shape can be found to accommodate octane. The 
pair functions goH(f) for bulk water and those around the hydration shell 
of water show that positions and widths of the two peaks are similar (see 
Fig. 4 in [39] . (3) The shape of the cavities in which the solute is housed are 
best described as oblate ellipsoids. Using A = j^-, where R3 and Ri are the 
major and minor radii of gyration of the solute, as a measure of anisotropy 
it was found that A varies from 1.0 for methane and ethane to 1.3 for octane 
[39] . (4) Several lines of evidence showed that octane and dodecane adopt 
extended structures in water. For collapse to occur a substantial fraction of 
dihedral angles must be in the gauche conformation. However, the probability 
of finding conformations with appreciable number of gauche conformations 
is small. In addition, the mean end-to-end distance is similar to that found in 
the gas phase, which implies that water does not force collapse of hydrocarbon 
chains even with n = 10. 

Collapse and association of Hydrocarbon chains in water A Flory- 
like theory showed that when n is ~ 20 a single HC would collapse and adopt 
compact structures [ID]. The estimate was made using average interaction 
per contact between the interaction centers, and hence is likely to be only 
approximately correct. Explicit simulations, first carried out by Mountain 
and Thirumalai [10J, confirmed that collapsed states can be realized with 
n = 18. However, subsequent studies [2,9,20,19] showed that, depending 
on the force field, realization of stable collapsed structures might require 
n > 20. Somewhat surprisingly, several studies have shown the HC in the 
collapsed structures are not merely maximally compact but adopt regular 
structures ranging from helix hairpin (HH), toroids and other complex folds 
depending on n [4T)ll9"ll2"0"] . Thus, linear hydrocarbon chains, which typically 
become compact on nanosecond time scales, are ideal model systems for 
probing some aspects of protein collapse. 
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We first report preliminary molecular dynamics simulations of the dynam- 
ics of collapse of a single hydrocarbon chain (HC) CH^ — (Cifa) 18 — CH3 
in water using the models described in the Methods section. To illustrate 
the pathways explored by the HC we show in Fig. [5^ the structures ex- 
plored in the transition from extended structures to a compact conforma- 
tion in a trajectory with first passage time on the order of 0.25 ns. Al- 
though, Fig. [2^ shows a sequential transition the movie (can be viewed 
at |http:/ /www.youtube.com/watch?v= dWW83YVhX^U| ) shows that tran- 
sitions between extended and collapsed structures occur prior to formation 
of the compact conformation. It appears that for collapse of extended HCs 
water expulsion occurs on time scales that are significantly less than r c , the 
overall collapse time. The radius of gyration of the final conformation of HC 
chain is consistent with Flory's scaling law ~ N u , where N is number of 
beads in the HC chain and v = 1/3, and with a pre- factor that is close to 
the LJ a. The structure of the compact structure is best characterized as a 
1.5 turn HH. 

In order to obtain a picture of HC association we also conducted prelimi- 
nary studies of interaction between two identical eicosane HCs. In the gener- 
ated trajectory we observed collapse of each HC chain and subsequent asso- 
ciation (see|http://www.youtube.com/watch?v=KP8bP6jWV40 for a movie 
of the trajectory). We should emphasize that this is only one of the possible 
pathways and more extended simulations are needed to fully determine the 
mechanism by which HCs associate. Similar to the case of single HC chain, 
the collapsed hairpin conformation can become extended (see the cyan chain 
in Fig. HJd), and re-collapses to the more condensed conformation, where the 
two HC chain tends to stack together with the long axis of one hairpin be- 
ing perpendicular to the long axis of the other. This resulted in a compact 
structure in which each HC adopts a HH conformation. It is likely that in 
this pathway water is expelled early in the association process, just as in 
oligomer formation between A/3 fragment with predominantly hydrophobic 
residues [2§], leading to the formation of a "dry" hydrophobic interface. 

3.2 Effects of Nanoconfinement: 

Water structure in nanopores: The water structure in the nanopores is 
considerably perturbed by the presence of nonpolar walls. Water hydrogen 
bonds are broken in the layer adjacent to the confining boundary result- 
ing in water being more loosely bound in the surface layer compared to the 
bulk phase. Water molecules are oriented with the molecular dipoles nearly 
orthogonal to the droplet radius leading to vapor-like, low density layer adja- 
cent to the walls, which as we describe below, can easily accommodate small 
hydrophobic species such as methane |60j . 

Methane is localized at the surface of a spherical cavity: Dis- 
ruption of water structure adjacent to the walls makes the surface region 
strongly favorable for nonpolar solutes like methane. Figs. [3] a and b show 
the free energy of a single methane molecule (black, filled circles) in spherical 
droplets of diameters 3 and 4 nm respectively. In both droplets the surface is 
more favorable for the methane than the interior by ~ 10 — 15 kJ/mol (4 — 6 
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k B T at 298 K) [gg pH H l BD ] . In the largest and the most bulk-like droplet, 
the free energy difference of ~ 10 kJ/mol is close to the hydration free energy 
of a methane molecule calculated from simulations and experiment |21| . This 
supports the view of the droplet surface inducing a vapor-like state of low 
density with broken hydrogen bonds [64] . Because of the broken H-bonds at 
the surface, the solvent loses less entropy due to ordering of hydrogen bonds 
around the methane if the nonpolar solute is confined to the droplet surface. 
The physical picture is supported by calculations, which show that the free 
energy of solvation of methane in spherical droplets is strongly dominated 
by entropic contributions [50] ■ Remarkably, our simulations and theory show 
that methane would be pinned at the surface regardless of the size of the 
spherical droplet. 

Ion solvation depends on sign of charge and size of the spherical 
droplet: The solvation free energies of positive (red, plus signs) and negative 
(blue, filled triangles) ions, M q ±, that are constructed by assigning charges 
of magnitude 0.4e to methane molecules arc also shown in Figs. |3J a and b. 
These fictitious ions have a greatly reduced surface propensity compared to 
methane. Unlike the nonpolar solutes, the free energies of the ions strongly 
depend on the droplet size. In droplets of diameter 3 nm and less, the ions 
are preferentially found at the surface. In the D = 4 nm droplet (Fig. [3Jd) 
the interior is nearly bulk-like [50] , and is more favorable than the surface for 
both positive and negative ions. Thus, the penetration of the ions into the 
droplet interior increases with decreasing curvature of the droplet surface, in 
accord with the findings for the chloride ion [53]. Just as in bulk water [31] 
there is an asymmetry between the behavior of cations and anions. Anions 
have a stronger tendency for surface solvation than cations with the same 
charge magnitude. As D increases, the enthalpy gain due to the solvation of 
the ions is greater than the entropy loss due to ordering of water molecules 
around this solute. The enthalpy-entropy balance depends on the interplay 
between the hydrophobicity and the ionic nature of the solutes. 

Hydration depends on charge density: To probe the balance between 
hydrophobic and electrostatic interactions, we calculated the free energies of 
ions, M q ±, with charge magnitude e, and also the ions M* q ± which have the 
same charge magnitude, but with twice the volume of a methane molecule 
(Fig. [3J:). Thus, each starred ion has half the value of v£ of the correspond- 
ing unstarred ion. The ions M q ± have a strong (enthalpic) preference for 
the interior over the surface of the droplet, with the free energy difference 
being of the order of 10 — 15 k B T . Molecular dynamics simulations with a 
non-polarizable force field [35] have obtained similar values for the the free 
energies of sodium and chloride ions in water, near hydrophobic or purely 
repulsive surfaces, and also at a water liquid-vapor interface. For M* q+ and 
M* g _, the surface is much less unfavorable compared to M q+ and M 9 _ re- 
spectively. The free energy profile for M* g _ illustrates the balance between 
hydrophobicity and ionic character. The anion M* g _ preferentially resides, 
not in the interior or at the surface, but 0.2 nm below the surface, due to its 
competing tendencies for surface and interior solvation. 

The asymmetry in the solvation of cations and anions in shperical water 
droplets, along with the £ dependence also reflects the differences in the sol- 
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vation of sodium and halide ions in water clusters and slabs |23| . Sodium and 
fluoride ions were found to be solvated in the interior while the larger halidcs 
had a propensity for surface solvation in the order Cl~ < Br - < I [53] • The 
surface propensities of solutes at an air-water interface depend on the bal- 
ance between hydrophobic (excluded volume) and charged interactions. This 
balance can be tuned by adjusting the ( (or by adjusting the van der Waals 
parameters, keeping the charges fixed). Decreasing ( increases the tendency 
for surface solvation. For a given £, anions have a greater preference for the 
surface than cations. Thus, a single parameter (£) tidily explains the nature 
of solvation of spherical ions at interfaces. 

Water-mediated interactions between solutes in spherical and 
cylindrical pores: Since methane molecules are confined to the surface of 
water droplets, their mutual solvent mediated interaction should reflect the 
disruption in the water structure due to the confining walls. Potentials of 
mean force (PMFs) between two methane molecules in droplets at various 
D values (Fig. [4^,) show an increased tendency for the methanes to associate 
compared to bulk water |49j . Surprisingly, we find that the secondary solvent 
separated minimum (SSM) that appears in the bulk profile is completely 
absent, even when D = 4 nm. The absence of the solvent separated minimum 
in the D = 4 nm droplet, which is entirely due to the presence of boundaries, 
is intriguing because the properties of water in large droplets arc bulk-like 
in all crucial respects |60| . The second minimum in the bulk corresponds to 
configurations where the methanes are separated by a single water molecule 
that is hydrogen bonded to other waters. In the immediate vicinity of the 
confining surface, where bonds between water molecules are disrupted, such 
configurations are strongly unfavorable. The preference for being pinned at 
the surface is consistent with a view that a methane is likely to be localized in 
regions with maximum unsatisfied hydrogen bonds, which is readily realized 
at the boundary. A striking feature of Fig. [4^, is that the calculated PMFs are 
independent of the size of the droplet. This is because methanes approach 
each other along the surface which is energetically similar in all the droplets 

Fig. HJd shows the PMFs between the ions M q ± with a charge magnitude 
of 0.4e. Reflecting the surface propensities, which depend on the droplet size 
(Figs. [5h. and b), the solvent- mediated interaction between the oppositely 
charged ions also varies strongly with the droplet size. As the size of the water 
droplet increases, the interior becomes more bulk-like and for the largest 
droplet, the profile is very similar to that for two methanes in bulk water, 
with a solvent separated minimum at a separation of ~ 0.7 nm. In the 2 and 
3 nm droplets, the secondary minimum is destabilized by <~ 2 kJ/mol relative 
to the 4 nm droplet, reflecting the greater role of confinement. 

The role played by the shape of the confining volume in determining the 
solvent-mediated interactions between two methane molecules is illustrated 
in Fig. [5j for a cylindrical water-filled pore of diameter D = 1.6 and length 
L = 2.08 nm. The calculated interaction free energy (solid, black line) shows 
the characteristic contact minimum and two distant minima corresponding to 
the pore dimensions |59j . Parenthetically, we emphasize that the free energy 
profile in the cylindrical pore cannot be directly compared to the PMFs in 
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bulk or in the spherical droplets (see Methods) . The distant minima arise due 
to the preferential solvation of the methanes at the surface of the cylindrical 
pore, just as in the spherical droplets. Treating the spherically symmetric 
methane molecules as point objects that are strictly conhncd to the pore 
surface, we also calculated the contribution from the translational entropy 
of the hydrophobes, —TAS A . This purely geometric term is shown by the 
dotted blue curve in Fig. [SJ As with the total free energy, —TAS A is only 
determined to within an additive constant, and is vertically aligned with the 
first distant minimum in the free energy at r w min{Z),L}. Fig. [5] shows 
that —TAS A qualitatively captures the curvature of the free energy profile. 
Thus, the solute translational entropy favors the first distant minimum in 
the free energy, while the solvent drives the solutes to cither minimize or 
maximize their separation. Qualitatively similar free energies of interaction 
have been calculated for the nonpolar pair of amino acid side chains ALA- 
PHE in cylindrical pores [STj . 

Balance between hydrophobic and electrostatic forces in LYS- 
GLU interactions: Many eukaryotic proteins are intrinsically disordered, 
which implies that they do not adopt well defined structures in isolation but 
do so only upon interaction with other partners. In general such proteins con- 
tain charges, and hence understanding water-mediated interactions between 
charged species both in bulk and in confined spaces determine their con- 
formational properties. The interactions between charged residues (Lys and 
Glu for example) also involve a subtle balance between hydrophobic and elec- 
trostatic interactions. The side chains of lysine and glutamatc, which carry 
charges of +e and -e, respectively, at physiological pH, also have nonpolar 
parts that prefer to be sequestered from the water. Therefore, the LYS + - 
GLU~ pair interaction can be mostly electrostatic or both electrostatic and 
hydrophobic, depending on their relative orientations. The PMFs between 
LYS and GLU in the —><!— and the tt orientations (sec Fig. [B^ and stick fig- 
ures in b and c) in the bulk show well defined contact and solvent separated 
minima. These configurations are stabilized in contact by 8.5 — 10 kJ/mol 
relative to extended separations. In the nanopores, the — pair, which is 
predominantly electrostatic, has the same contact and solvent separated min- 
ima at ~ 0.6 and ~ 0.8 nm respectively, as in bulk solvent (Fig.JBJ)). In this 
orientation, the free energy has another minimum corresponding to the state 
where the two molecules are approximately a pore length apart, with their 
charged ends hydrated and the nonpolar ends at the surface. The free energy 
in the tt orientation shows features of both electrostatic and hydrophobic 
interactions. This pair is strongly driven to contact at ~ 0.43 nm. 

The LYS- GLU system illustrates how confinement alters the balance be- 
tween hydrophobic and electrostatic interactions. In the absence of charges, 
the interaction between these side chains would be entirely hydrophobic. If 
this were the case, then the free energy of interaction between these hypo- 
thetical molecules would have two distant minima in addition to the contact 
minimum, at separations corresponding to the pore dimensions, similar to the 
case of the two methanes and the ALA-PHE pair [BT|. If this was a purely 
electrostatic pair, the solutes will be fully hydrated and the free energy pro- 
file will have a contact minimum and an SSM, but no distant minima. In this 
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scenario, the presence of the confining boundaries will have only a minimal 
effect on the interaction between the solutes. Because favorable electrostatic 
interactions are expected to dominate in the bulk, we expect that in most 
cases LYS and GLU will retain their hydration shells. In accord with this 
observation, it is found that LYS and GLU in A/3-peptides are solvated in 
bulk water }271 I551. The orientation dependent interactions in the nanoporc 
show that confinement (Fig. [SJa) can promote salt bridge formation between 
LYS and GLU. Fig. [St shows that the LYS-GLU system is largely insensitive 
to the solvent density. Even with the net solvent density reduced by 50%, 
both the — and the tt interactions are similar to those at bulk solvent 
density (Fig. [BJd). The distant minimum at ~ 2.3 nm for the — ><— orientation 
disappears with the decrease in solvent density. 

Confinement enhances peptide stability: Our primary motivation 
for undertaking these studies was to gain insights into the stability of pro- 
teins in confined spaces. The extent of confinement-induced stabilization (or 
dcstabilization) of peptides or proteins, with respect to bulk, will depend on 
the peptide sequence, entropy loss of the denatured state ensemble (DSE) 
and solvent density. (1) While nonpolar SCs prefer to be localized near the 
hydrophobic pore boundary, oppositely charged side chains can be stabilized 
in contact either by the formation of a salt bridge or by interior solvation. 
Confinement also drives polar SCs to form contacts. (2) Confinement in a 
narrow cylinder, as is the case in ribosomes, also restricts the rotamer degrees 
of freedom, thus selecting a particular orientation that may not be preferred 
in the bulk. The decrease in the rotamer degrees of freedom in nanopores, as 
well as overall restrictions in the conformations of the peptide increases the 
entropy of the DSE. (3) Confinement destabilizes the contact minimum be- 
tween hydrophobic SCs only if the water density in the nanoporc and the bulk 
are similar. At lower densities the contact minimum is favored |61j . These 
observations suggest that for a generic polypeptide sequence, confinement 
should enhance the stability of ordered states. 

To test these predictions, we simulated two tripeptides, (Ala) 3 , a hy- 
drophobic sequence and Lys-Ala-Glu, a hydrophilic sequence with no net 
charge, in bulk water and in two cylindrical pores (Fig. [7]). Both peptides are 
in their zwitterionic form, terminated with a positively charged amino group 
and a negatively charged carboxylic group. In the pores, we find that the 
center of mass of the (Ala)3 backbone is close to the surface. The charged 
backbone is hydratcd, while the nonpolar side chains are sequestered from 
the water. Not surprisingly, the center of mass of the hydrophilic sequence 
remains close to the pore axis, away from the surface. These findings are in 
accord with our expectations based on interactions between amino acid side 
chains in cylindrical nanopores. 

The changes in stability of the peptide conformations, with respect to 
the bulk, can be assessed by computing the distributions of the backbone 
dihedral angles <j> and tp of the middle residue in each peptide (Fig. HJ). The 
a-helical region is defined as -80° < 4> < -48° and -59° < tp < -27° 
[27] . We find that cylindrical confinement stabilizes both peptides, with the 
hydrophilic sequence being more strongly stabilized (See Table[l|. The extent 
of stability clearly depends on the sequence and its length, pore geometry 
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and solvent density, just as predicted using the results for PMFs between 
amino acid side chains. 



4 Conclusions 

Model systems, such as hydrocarbon chains that can adopt compact struc- 
tures in water, provide useful caricatures of the earliest stages of collapse 
transition in proteins. It is now firmly established that only for hydrocarbon 
chains longer than n ~ 20 the collapsed structure is stable. The preference for 
the bonds to be in all trans configuration prevents smaller hydrocarbon chains 
(at least sizes on the order of dodecane) from adopting globular conforma- 
tion. Surprisingly, the collapsed structure is ordered. The chains adopt HH 
and higher order structures such as toroids depending on the length of the 
hydrocarbon chains. Similarly, association between two hydrocarbon chains 
produces an ordered structures in which each chain adopts a 1.5 turn HH. 
It appears that in both processes water expulsion occurs early (before the 
chains adopt regular structures) although much more exhaustive simulations 
are required to quantify the time scales and the mechanisms of assembly. 

Ion solvation in confined space not only reveals an asymmetry between 
positive and negative but also shows the importance of charge density, £, 
in the preference for interior and surface solvation. Charge density, a single 
parameter, accounts for excluded volume and electrostatic effects. Our work 
predicts that lower charge density ions prefer the surface and as the charge 
density increases interior solvation occurs. For a given £ anions have greater 
preference for the surface than cations. Most studies (sec for example [23] ) 
have argued that ion polarization has to be accounted for in order to explain 
the enhanced propensity of I -1 to be at the air/ water interface than CI -1 . 
Although polarization and size are not unrelated our work qualitatively ex- 
plains in terms of £ alone the tendencies of ions of differing sizes to localize 
near the surface. 

An important aspect of the current perspective is to provide molecular ba- 
sis for explaining stability of encapsulated peptides in cylinders. First, side 
chain orientation can dramatically alter interactions in confinement com- 
pared to bulk. Second, interactions between charged residues (lysine and 
glutamic acid for example) in confined space depends strongly on a balance 
between hydrophobic and ionic interactions. Using these results and loss of 
greater conformational entropy of the unfolded states relative to the folded 
structures upon confinement we have shown that for a generic amphiphilic 
sequence folded states should be stabilized when trapped in cylindrical pores. 
Explicit simulations in water |61j and extensive coarse-grained simulations 
of peptides in carbon nanotubes [U| support the theoretical expectation. 
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(a) 




Fig. 1 (a) Stick representation of a HC chain (CH 3 - {CH2) 18 - CH 3 ) (b) Each 
methane group is represented with a bead centered at the carbon atom. 




Fig. 2 (a) Representative snapshots of a trajectory showing a single HC chain 
collapses in water (b) Snapshots of a trajectory showing two HC chains collapse 
and associate in water. In both cases the HC chains adopt ordered helical hairpin 
structures. 

Table 1 Stability changes in peptides confined to cylindrical pores 



Pore size (D, L)^_ 


AF [k s T^ 


AF [k s T^ 


2.0, 2.9 


-0.61 


-0.84 


1.4, 2.9 


-0.14 (-0.19^) 


-1.50 (-1.02 4 ) 



a D and L axe the diameter and length of cylindrical pores in mn. 
b AF = — log (P a /Pa) f° r (Ala) 3. P a and are the probabilities of being in 
the a-helical region in the pore and bulk respectively. 
c AF for Lys-Ala-Glu. 
d AF for p = 0.5p bulk 
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Fig. 3 (a,b) Free energies of methane molecules and ions M 9 ±, derived by as- 
signing charges of ±0.4e to the methanes, in spherical water droplets of diameter 
D nm. The zero of the free energy scale for every droplet is at the surface. The 
methanes have a strong preference for the droplet surface at all values of D. In the 
case of the ions, the surface propensity depends strongly on droplet size. In every 
droplet, anions have a greater preference for the surface than cations of the same 
charge magnitude, (c) Free energies profiles for ions M 9 +, M g _, M* 9 + and M* 9 _ of 
charge magnitude l.Oe. The starred ions have twice the volume and therefore half 
the charge density of the corresponding unstarred ion. Curves for M 9 +, M 9 _ and 
M* 9 + are referenced to the origin (droplet center) while that for M%_ is referenced 
to the minimum at 1.3 nm for clarity. In all cases, only differences within each curve 
are relevant. 
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Fig. 4 (a) Potentials of mean force (PMFs) between two methane molecules in 
spherical droplets of diameter D nm. Just as in the bulk, there is a distinct primary 
minimum. However, the characteristic solvent separated minimum is absent even 
at D = 4 nm. The curves are shifted vertically so that the zero of the free energy 
scale is at contact for the two methanes, (b) PMFs between M q + and M q _, with a 
charge magnitude of 0.4e, in droplets of different sizes. 





5 




-1 

0.5 1 1.5 2 2.5 

r [nm] 



Fig. 5 Interaction free energies of two methane molecules in a cylindrical nanopore 
of diameter D and length L nm, containing water at bulk density. The dotted 
blue line shows the contribution from the translational entropy of the hydrophobic 
species, assuming that they are strictly confined to the pore surface. 
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Fig. 6 (a) Potentials of mean force between LYS and GLU in bulk water in the 
— ><— (stick figure in b), and tt (stick figure in c) orientations. Curves are translated 
vertically so that the zero of the free energy scale is at contact. The snapshot shows 
the — pair at the solvent separated minimum, (b) Free energies of interaction of 
LYS and GLU in a nanopore of diameter D and length L. Water is at bulk density. 
The inset shows the — >•<— pair at the distant minimum at ~ 2.3 nm. (c) LYS-GLU 
interaction free energies in the same nanopore at half the bulk water density. The 
inset shows the tt pair a t contact. 
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